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Magnetic oscillations and field induced spin density waves in (TMTSF) 2 C10 4 
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We report an analysis of the effects of magnetic field on a quasi-one-dimensional band of inter- 
acting electrons with a transverse dimerizing potential. One-particle problem in bond-antibond 
representation is solved exactly. The resulting propagator is used to calculate the spin-density-wave 
(SDW) response of the interacting system within the matrix RPA for the SDW susceptibility. We 
predict the magnetic field induced transition of the first order between interband SDWo and intra- 
band SDW± phases. We reproduce the rapid oscillations with a period of 260 Tesla and the overal 
profile of the TMTSF 2 C10 4 phase diagram. 
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Investigations of quasi-one-dimensional electronic sys- 
tems at high magnetic fields and at low temperatures con- 
tinue to give an important insight into the one-particle 
properties and interaction-induced phases such as spin- 
and charge-density-wave, superconductivity, and Mott 
localization. One of most spectacular phases of this 
kind are field-induced spin density wave (FISDW) , found 
in Bechgaard salts p| and in some other low-dimensional 
compounds. The phenomenon of the FISDW is well 
understood in the Bechgaard salt (TMTSF) 2 PF 6 where 
the cascade of SDW phases with quantized wave- vector is 
induced by orbital effects of magnetic field to the quasi- 
one-dimensional orbits of band electrons. Theory based 
on the mechanism of quantized nesting [Ej reproduces 
satisfactorily main experimental data for this salt. 

In this letter we concentrate on (TMTSF) 2 C10 4 , a 
Bechgaard salt which after a slow cooling || [| enters 
into a qualitatively different type of FISDW phase at 
low temperatures, with a phase diagram that is still, after 
more than ten years of intensive studies |^| , a matter of 
both experimental and theoretic controversies. In partic- 
ular for magnetic field B>8T the nature of the ordering 
in the relaxed material is not a simple FISDW with some 
low integer quantum number N, but a qualitatively difer- 
ent state containing several puzzling subphases.fg, |7[ ||] 
This phase is at 8T separated by a line of first order 
transition from a cascade of FISDW phases which very 
much resembles to that of the standard model. Another 
characteristic phenomenon, the rapid oscillations (RO) 
in 1/B with a frequency of 260 Tesla are visible in trans- 
port properties in both metallic and FISDW state. |^|, ||, |j| 
Similar RO are seen also in thermodynamic quantities 
like torque, magnetization, sound velocity and specific 
heat, but only in the ordered phase. [|[ || The highest 
value of T c in the T C (B) dependence is 5.5K, instead of 
12K as expected from analogy with the (TMTSF) 2 PF 6 
salt. 

The incompatibility of above facts with the quan- 



tum nesting model (QNM) for a single quasi-lD band 
is believed to stem from the particular ordering of CIO4 
anions. || This ordering introduces the new modulation 
with the wave vector (0,7r/6, 0), i. e. a dimcrization 
in the low-conducting direction with the inter-chain dis- 
tance b. The magnitude of the dimerizing potential can 
be tuned to some extent by varying the cooling rate 
|^|, |6). Thus, anions presumably remain disordered in the 
rapidly quenched samples. Then there is no dimerization 
gap in the band, and the system shows properties of a 
single quasi-lD imperfectly nested band with a SDW or- 
der appearing already in the zero magnetic field [B @, [[(J . 
The anion ordering in slowly relaxed samples is at about 
24K, and coincides with the onset of rapid oscillations in 
the magnetoresistance. |n| 

The dimerized band has two pairs of Fermi sheets in 
the new Brillouin zone. Already simple geometric ar- 
guments |?J suggest three possible nesting wave vectors 
favoring various SDW phases. First, interband nesting, 
leads to SDWo that is the two-band version of the stan- 
dard FISDW phase. Other two nesting vectors relate 
Fermi sheets within the same band. They give SDW+ 
for antibond nesting and SDW_ for bond nesting. How- 
ever the interplay between SDWo and SDW± is not only 
a geometric question of the choice of the nesting vec- 
tor. Due to a finite anion potential V in the kinetic part 
of the Hamiltonian an off-diagonal term appears in the 
SDW response, making necessary an appropriate matrix 
approach jl^, |l3) in the calculation of the critical suscep- 
tibilities. The response matrix is formulated in the space 
of two order parameters, ("homogeneous") and A a 
( "alternating" ) , determining the magnetic pattern 



in 



: (x, R±) = {A h ± A a ) cos [{2k F + k)x + pnd] . (1) 



Here d = 2b and the upper and lower sign stay for even 
{R± = rid) and odd {R± = nd + d/2) chains respectively. 
As it is shown in Refs. @ [T§, SDW [A h ^ 0, A„ = 0] 
is stabilized for low values of V (providing the imperfect 
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nesting parameter t' b allows for SDW stabilization) , while 
SDW± with A h > A a ^ gets stable for V/t b > 1.6 ir- 
respectively to the value of t' b . Here 4 is the interchain 
hopping integral. The slowly relaxed (TMTSF) 2 C10 4 
samples are expected to lie in the range of intermedi- 
ate values of V in which there is no SDW ordering at 
B = down to T = 0. 

Indeed, as it will be shown, V/tb fitting the experi- 
ments is close to one. In this range it is not allowed 
to use the quasi-classical approximation of Gor'kov and 
Lebed ]lq| , which consists in making Peierls substitution 
p — > p— eA in each sub-band separately and including the 
anions' effects only via magnetic breakdown (MB) junc- 
tions near the zone boundary. While this approximation 
is sufficient for V/tb <C 1 , here one has to solve the whole 
quantum-mechanical problem instead. 

It was pointed out several times (see Q and references 
therein) that a mechanism of coherent inter-band tun- 
nelling, very similar to Stark over-gap quantum interfer- 
ence (QI) in magnesium [Q, is essential for high-field 
physics in (TMTSF) 2 C10 4 . In particular, RO in metallic 
state can be explained only in terms of QI mechanism be- 
cause no closed orbits exist. On the contrary, in the SDW 
state both closed orbits and Stark interference contribute 
to RO. Oscillating behavior periodic in 1/B can be seen 
already at the level of one-particle spectrum. Namely, 
electronic propagator with longitudinal momentum k has 
poles at 



E f = v f [f(k - k F ) + GN] ± v F G5 , 



(2) 



where / is left-right index, G = eBb/h is the magnetic 
wavenumber and N is integer number. The first term in 
eq. (|^) is the standard QNM dispersion and the last term 
is the splitting due to anions. Over-gap resonances are 
present in 5(B) as will be discussed below (see Fig.0). 
The expression for the spectrum (|^) is common to per- 
turbation calculations [la], to quasi-classical tunnelling 
analysis |ll|, and to our exact solution as well. What 
changes from one approach to another are the depen- 
dence 5(B) and the result for electronic wavefunction. 
In order to obtain them exactly we start from the ef- 
fective one-particle Hamiltonian for electronic operators 
Vf(x,p) 

H Q = iv F p 3 d x + r 3 T(pb - Gx) + f(pb - Gx) - V n , 

(3) 

where p's and r's are Pauli matrices in left-right and 
bond-antibond indices respectively. The most general 
transverse dispersion was split into two parts 

oo oo 

T(pb) =2j2h cos[(2j - l)pb] , f(pb) = 2 ^ cos[2jpb] 
j=i j=i 

(4) 

corresponding to effective hoppings Jl7[ between odd and 
even neighbors respectively. We diagonalize Hq by the 



unitary transform 



9 4 



a f /3 f 



oif»if, f 



(5) 



with \a\ 2 + \{3\ 2 — 1, and functions a, (3 and 8 depending 
on x and p only through the combination z — pb — Gx. 
From the requirement that the effective Hamiltonian for 
field $ is only ifv F d x we get 



0(z) 



1 

V F 



dzf(z) 



(6) 



and a system of differential equations for functions a and 



ifv F a' f (z) = -T(z)a f (z) - Vf3* f (z) 
ifv F p' f (z) = -T(z)P f (z) + Va)(z). 



(7) 



Note that 8(z + 2tt) = 9(z) and that a + (z) = a*_(z) 
and (3+(z) — j3*_(z), so that it suffices to follow e. g. 
solutions a + (z), P + (z) of the system (^). According to 
Floquet theory these solutions can be written in the form 
a(z) = A(z) exp(— izS); (3(z) = B(z) exp(iz5). A and B 
are periodic with the period 2ir, and the closer inspection 
shows that the Floquet exponent 5 for the system (Q) is 
real for all values of parameters, at least after keeping in 
T(z) only the leading term t\. 

Once we find A, B and 5 the wave functions (x,p\Fk) 
of the states created by ^(x,p) are known. The cor- 
responding spectrum is one-dimensional, Ek — v F f(k — 
k F ). Projection of \Fk) to the plane wave (k',p\ is 

Y L k e^ N -^ pb 5[k' — k + G(N - 5)} ( °f 
+ Vk e* N+ V pb 5[k' - k + G(N + 5)] ( b & N ] | (8) 

Coefficients apf, bpf, &n and 6jv are Fourier compo- 
nents of the products Aexp(i8), Bexp(i9), A* exp(i8) 
and B* exp(i9) respectively. Coefficients fik and Vk are 
fixed by initial conditions. The expression (|§|) tells us 
how the plane wave exp(ikx + pR±) is decomposed into 
discrete states N. Each state N is split by 5 in a way 
that components with the tilt —5 have the statistical 
weight |^.| 2 (|ajv| 



+ \b N \< 

12/ 



and the ones with the tilt 
-5 have the weight | ^ fc I ^ ( I ^ at | 2 + |aAr| 2 ). Green function 
(&(x, p)^ (x' , p)) is easily constructed using transforma- 
tion (^) and knowing that ($$^) = (iu n — iv F fd x )^ x = 

The Floquet exponent 5 and the functions A and B 
are calculated using the Hill's theory and the fundamen- 
tal matrix method Jll| . In the present work we limit 
our calculations to first harmonics in Eq. |^ only, param- 
eterized with t\ = tb and t[ = t' b . Let us concentrate on 
the magnetic field dependence of the Floquet exponent <5 
that splits the QNM spectrum as given by Eq. S. Figure 
|l|(a) shows the energy ui c 5 (in units of V) as a function of 
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the magnetic breakdown parameter re = 2uj c tb/V 2 , where 
u c = vfG is the cyclotron frequency. In quasi-classical 
picture re determines the probability of the over-gap tun- 
nelling P — exp(— 7r/2/c).E^| One sees that the crossover 
from oscillating to saturating behavior does not coincide 
with the crossover from the weak (re < 1) to the strong 
(re > 1) MB. The position of the last zero of 8 is not uni- 
versal in re, but approximately in r = \(^V) 2 + tl] 1 ^ 2 /uj c , 
where the value of 7 is 0.77. Fig. [j](b) shows S(r) for 
several "polar angles" defined by tan# = tb/-fV '. Os- 
cillations of (5 are approximately periodic in r with a 
period of 0.80. Choosing the parameters tb — 300-ff, 
v F = 2 x 10 5 m/s, and b = 7.7x 10" 10 m we fit RO at 260 
Tesla by putting V « 0.8t b . 

Taking the limit of strong magnetic filed uj c /tb 3> 1 and 
of weak anion potential V/tb <C 1 we can easily reproduce 
the ID spectrum of Osada et al. p5|, —> fvp(k—kf)± 
lo c 5 with 8 — > (V/Lo^J^Atb/Ldc), Jo being the Bessel 
function. On the other hand the spectrum of Gor'kov 
and Lebed [|16[ is reproduced for weak anion potential, 
V/tb <C 1. [|18|| The above fit, as well as other insights |2(], 
|l| however strongly suggest that V in (TMTSF) 2 C10 4 
is rather large, i. e. comparable to 

We proceed with the solution of the interacting prob- 
lem. Neglecting the absence of a presumably small Umk- 
lapp scattering, the effective coupling for SDW is the 



forward scattering amplitude 32, here simply denoted by 
U. We employ the matrix RPA formalism developed 
in Ref.|12|. The resulting relevant bare susceptibility is 



Xl(q; T) = \{Xaa + Xhh + [(Xaa - Xhh) + MXha) 2 ] 1/2 }, 

playing in the Stoner criterion 



l-U X i(cLc,T c )=Q, 



(9) 



q c being the wave vector at which Xi(q) has the max- 
imum. The ratio of two SDW order parameters from 
Eq.(|l]) is also a function of bare correlators Xaa, Xhh, Xah 
in the (a,h) basis (see |l2|). The bare correlators in the 
magnetic field are given by 
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(10) 



where P ,P± stand for P(q\\ —NG,T) and P[q {] -G(N± 
28), T] respectively, P(k,T) being the familiar ID Lind- 
hard function at the wave number 2k F + k. Pq and P± 
are the inter-band and the intra-band susceptibilities of 
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the N-th split level in the decomposition (||). The de- 
pendence on the transverse momentum is present in the 
amplitudes I(q±,N), 

I h0 (q±,N) = J2 (a»b N - n - 6 n ajv-„) e l( - n ~ N / 2 ^ 

n 

I h+ (q±,N) = J2 {Ka N -n + Mtf-n) e ^ n - N ^ 

n 

h-(q±,N) = J2 (anajv_„ + b n b N ^ e ^ n ~ N ' 2 ^ 

n 

I a0 (qx,N) = (a n a N - n - b n b N ^ e i( - n ~ N ' 2 ^ 

n 

I a+ (qx,N) = {a n b N -n + MiV-n) e 1 ^' 2 ^ 

Th 

n 

(ii) 

There are two important selection rules for these ampli- 
tudes, namely for N even, Iho(N) = I a o{N) = while 
for N odd, I h ±(N) = I a ±(N) = 0. Thus the interband 
processes contribute only to FISDW phases with odd N 



while the intraband processes contribute only to phases 
with even N. Consequently only phases with even N 
"see" the splitting by 5. 

The q— dependence of the susceptibility Xi(q) for 
the particular choice of parameters, uj c = O.ltb.V — 
0.8t b ,t' b /t b = 0.03, T/t b = 0.001, is shown in Fig. |. 
The overall envelope assumes the shape present already 
in the absence of magnetic field |L2|. It is now superim- 
posed by a well known characteristic of FISDW suscep- 
tibilities |U, logarithmic peaks corresponding to single 
one-dimensional bubbles P(k), weighted by p-dependent 
amplitudes as defined by Eqs. (|l0| , |ll| ). Qualitatively new 
feature regarding these peaks is the splitting of peaks 
with even N by ±5 around the positions at k = NG. 

According to Eq.(^) at T — T c the highest of peaks 
in Fig. ^ attains the value 1/U. Fig.|| shows the result- 
ing phase diagram for a realistic choice of parameters, 
V = 0.85*6, t' b = 0.03t b and T C (V = t' b = 0) = 13K. The 
resulting maximal critical temperature within the present 
field range is T c ! Iiax ps 1.1 K. The most obvious charac- 
teristic of the obtained phase diagram is the first order 
transition from SDWo to SDW± at B c ~ 9 Tesla. Depen- 
dence T C (B) for B < B c is similar to the FISDW cascade 
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FIG. 3: (a) Phase diagram, (b) Energy ratio u> c 5/tb on the same magnetic scale as the phase diagram. 



in TMTSF 2 PF 6 , with the difference that here only odd 
phases appear because the even ones are suppressed by 
splitting. For B > B c the critical temperature increases 
towards the highest value T c max . As the magnetic field 
further increases the critical temperature T C {B) starts to 
oscillate, with the sharp dips corresponding to commen- 
surability condition 2GS = G between the Floquet wave 
number and the magnetic wave number. 

Note that the present choice of anion potential V — 
0.85t b places us in the intermediate regime on scale V 
[HU where the response of both SDWo and SDW± is sup- 
pressed for B = 0. Hence, being of comparable magni- 
tude the two instabilities compete once they are restored 
by magnetic field. Both phases in fig.||are sensitive to V , 
general trend being that by increasing V one decreases 
T C {SDW ) and increases T C (SDW±). The parameter of 
imperfect nesting in the standard model, t' b , here affects 
only SDWo , while for SDW± t' b plays a role of an effective 
nearest neighbor hopping. We remind that the effective 
parameter of imperfect nesting for SDW± is a function 
of t b /V, as pointed out in Ref. jL2l. t' b acts on SDWo in 
a standard way 0, i. e. it fixes the width of the FISDW 
cascade Au> c ~ t' b so that by increasing t' b one reduces 
T C (SDW ) at fixed B. 

The result of the subtle interplay between two scales 
V and t' b is that the realistic profile of the phase diagram 



is possible only within a rather restricted range of the 
(V,t' b ) space. By increasing V or t' b by a few percent 
one reduces T C (SDW ) below T C (SDW±) in the whole B 
domain. On the other hand by decreasing V by a few 
percents one gets a hump in T c (SDWq) on the left of the 
transition SDW - SDW±. 

The rapid oscillations in observable response func- 
tions are related to the oscillations of 5 shown in 
Fig.^(b). Generally, we expect RO to be visible if two 
conditions are fulfilled. First, in order to have an over- 
gap interference one needs a moderate MB parameter, 
k ~ 1. Second, one has to be in the oscillating regime 
on the scale r, which is equivalent to n < p(V/t b ) = 
2(V/t b )- 2 yJ(^V/t b ) 2 + 1. FigJ|(b) shows how the energy 
tilt u> c 5 varies in various parts of the phase diagram. At 
30 Tesla we have n ~ 0.5 and p(0.85) w 3.3, so that both 
conditions are fulfilled. The effect is expected to be even 
stronger for higher fields because k then increases. 

The maximal value of the critical temperature in FigJ^, 
ymax K i j s considerably smaller than the experi- 
mental value of 5.5K. In this respect we note that T c max 
is essentially model dependent quantity, i. e. that the 
Hamiltonian (||) represents a minimal model for under- 
standing the interplay between two SDW phases in the 
magnetic field. Namely, recent experiments sug- 
gest that the anion ordering in TMTSF2CIO4 induces 
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also, beside a strong dimerizing potential V, rather large 
changes in other band parameters. 

The present treatment also does not include the quan- 
titative analysis of the splitting of degeneracy of two 
intraband phases, SDW+ and SDW-. Physically the 
degeneracy is lifted because the realistic tight-binding 
dispersion along the chain is not strictly linear. Conse- 
quently the dominant instability will be that of SDW—, 
as discussed in Ref. jl^|. Similar conclusions were ob- 
tained also by numerical calculations p2| ], but without 
taking into account the two component aspect of the or- 
der parameter The second critical temperature can 
be calculated within Landau theory as in Ref.|ll|, and 
by taking the nonlinearity of the band dispersion into 



account. The subphases of the high-field phase corre- 
spond to SDW+ phases within SDW-, each one nesting 
its own pair of Fermi sheets. Such scenario is impossible 
for SDWo since it proceeds through nesting of all four 
sheets at the single critical temperature. 

In conclusion, we have solved exactly the one-particle 
problem of dimerized Q1D band of electrons in magnetic 
field. Observables contain characteristic periodicity in 
1/.B, consistent with 260 Tesla oscillations in normal and 
SDW phases of (TMTSF) 2 C10 4 . Using matrix RPA for 
SDW susceptibility we reproduce the first-order transi- 
tion between two types of FISDW ordering, as well as 
the overall profile of the experimental phase diagram. 



[1] 
[2] 



[3] 
[4] 



[5] 
[6] 

IT] 
[8] 
[9] 
[10] 

[11] 
[12] 

[13] 



T. Ishiguro, K. Yamaji, and G. Saito, Organic Supercon- 
ductors He, (Springer- Verlag, Berlin, 1998). 
P. M. Chaikin, J. Phys. I (France) 6, 1875 (1996); P. 
Lederer ibid. 6, 1899 (1996); V. M. Yakovenko and H. S. 
Goan ibid. 6, 1917 (1996). 

N. Biskup et al, Phys. Rev. B 60, 15005 (1999). 
For historical references and recent developments of the 
quantum nesting model, see A. G. Lebed, Phys. Rev. 
Lett. 88, 177001 (2002). 

J. S. Quails et al, Phys. Rev. B 62, 12680 (2000) 

N. Matsunaga et al, Phys. Rev. B 62, 8611 (2000); Icond 



mat/0206010 



S. K. McKernan et al, Phys. Rev. Lett. 75, 1630 (1995). 
O.-H. Chung et al, Phys. Rev. B 61, 11 649 (2000). 
W. Kang et al, Synthetic Metals 120, 1073 (2001). 
A. Bjelis and K. Maki, Phys. Rev. B 45, 12887 (1992). 
S. Uji et al., Solid State Commun.103, 387 (1997). 
D. Zanchi and A. Bjelis, Europhys. Lett. 56, 596 (2001). 
K. Sengupta and N. Dupuis, Phys. Rev. B 65, 035108 



(2002). 

[14] R. W. Stark and C. B. Friedberg, J. Low Temp. Phys. 

14, 111 (1974). 
[15] T. Osada et al, Phys. Rev. Lett. 69, 1117 (1992). 
[16] L. P. Gor'kov and A. G. Lebed, Phys. Rev. B 51, 3285 

(1995); ibid. 1362. 
[17] K. Yamaji, J. Phys. Soc. Jpn. 51, 2787 (1982). 
[18] D. Radic, A. Bjelis and D. Zanchi, to be published. 
[19] S. L. Ross, Diferential Equations, 3rd edition (John 

Wiley & Sons, New York, 1984). pp. 505-521; E. 

L. Ince, Ordinary Differential Equations (Dover Publ., 

1956) pp. 384,503,507; E. Kamke, Differentialgleichun- 

gen (Akademische Verlag. Becker & Erler Kom.-Ges., 

Leipzig, 1942). 
[20] H. Yoshino et al, J. Phys. Soc. Jpn. 68, 3142 (1999). 
[21] D. Le Pevelen et al., Eur. Phys. J. B19, 363 (2001). 
[22] K. Kishigi, J. Phys. Soc. Jpn. 67, 3825 (1998). 



